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Abstract 

We study a class of models for neutrino mass matrix in Type-II seesaw with A 4 family symmetry. The 
resulting neutrino mass matrix can be naturally made to respect a // — r exchange plus CP conjugate 
symmetry (GLS) with the CP violating phase 6 and the mixing angle 023 predicted to be ±7r/2 and vr/d, 
respectively. When GLS is explicitly broken by complex Yukawa couplings, the model predictions for 6 
and 023 can be significantly modified. Should future experiments will indeed determine 623 and 6 cp away 
from the GLS limit values, one then had to consider models with broken GLS. We study several simple 
scenarios to show how the modifications arise when GLS is broken and how future experiments can test 
this class of models. 
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I. INTRODUCTION 


Experiments have collected a lot of precious information about the neutrino mixing parameters. 
The mixing angles in the Pontecorvo-Maki-Nakagawa-Sakata matrix[l| Vpmns are not always small 
as their quark mixing counter partj^, |^. In the standard parametrizationj^, for three 
neutrino mixing commonly usedj^, |^, the mixing angle 623 is close to 7r/4, 612 is also large, ^13 is 
relatively small but away from zero. There are also evidences showing that the CP violating Dirac 
phase 6 is close to —7r/2 (or 37r/2). If these data are further conhrmed, the neutrino mass matrix 
will have a simple form. Assuming that neutrinos are Majorana particles, the neutrino mass matrix 
dehned by the term giving neutrino masses in the Lagrangian (l/ 2 )himj,n£ has the following form, 

= ^PMNS'^y^PMNS i (^) 


where rh^ = diag{mi,m2,m3) with rrii = \mi\exp{iai). Here we have put Majorana phase informa¬ 
tion in the neutrino masses. The standard form for Vpmns is given by 


VpMNS — U (6^12, 6*13, 6*23) 


I '^12 ^23 
\ '5i2'523 


^12^13 

-C 12 S 23 S 13 

“ '^12'^23'® 13 


S 12 C 13 SigC 

^12 ^23 ~ '5i2 ^23 '^13 '^23 ^13 

“*^12 ^23 ~ '®12 ^23 *^23 


( 2 ) 


where Cij and Sjj are cos6*jj and sin^jj, respectively. They are all normalized to be positive. 
With 5 = —7r/2 and 6*23 = vr/d, rrii, has the following form[ 3 , 0 , [lO 


rriy 


a c + ijS —{c — i( 3 ) 

c + i /3 d + i'j b 

— {c — i/ 3 ) b d — 


( 3 ) 


where 


a = miCi 2 Ci 3 m 2 S? 2 C ?3 - 7713^13 ) 


C =- - m2)si2Cl2Ci3 , 

v2 

/3 = -^Si3Ci3 (miCi2 m2Si2 + m3) 


^ (’^l(s?2 + C?2S?3) + ’^2 (c?2 + SI 2 SI 3 ) - m3C?3) , 

^ ~ 2 ~ C 12 S 13 ) -|- 772 - 2(012 — ■S 12 S 13 ) -|- 7723013 ) , 

7 = -(7721 - 7722 )Si 20 i 2 Si 3 . ( 4 ) 


Note that in the most general case, because non-zero Majorana phases, the parameters a, b, o, d, 
/3 and 7 are all complex. 

One has the degrees of freedom to redehne the neutrino helds phases and the most general form 
of the above mass matrix can be rewritten as 

0 0 \ / a c + i (3 -{c-i/ 3 ) \ / 0 0 \ 

rriiy = I 0 0*^2 ^)( c-\- i /3 d + i'^ b )(^ ^ I ’ 

\ 0 0 y I d-z^ J \ 0 0 y 

where the phases pi are arbitrary. One can choose some particular values for pi to obtain forms 
of 772 j, for convenience of analysis. For example the sign for the “ 13 ” and “ 31 ” entries can be 
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removed by choosing pi = P 2 = 0 and = vr, the resultant matrix can be written in a more familiar 
forms 


/ a c + ip {c — iP) \ 

= \ c + ip d + i'^ b , (6) 

y (c — iP) b d — i'y J 

where b = —b. 

If neutrinos do not have non-trivial Majorana phases (but mass can be positive or negative), 
the mass matrix has the following form 

/A C C*\ 

m, = \ C D* B \ . (7) 

\C* B D ) 


Replacing 5 = — 7 r /2 by 5 = 7 r/ 2 , the neutrino mass matrix is given in a similar form as that 
in eq-dH]), but P and 7 need to be multiplied by a sign. This implies that without further 

information given, a general mass matrix in the form given by eq.(l 6 ]) can give 6 = ±7r/2 and 
^23 = 7 r/ 4 . Whether they predict + 7 r /2 or — 7 r/ 2 , additional information need to be provided [lll|. 
It has been pointed out that the general form in eq.([n]) is a necessary condition for 6 = ■ii'K/2 and 
^23 = 7 r/ 4 , but not sufficient condition. In our later discussions, unless specihed, the mass matrix 
of the form in eq. ([H]) is always referred to a general form whose elements are not necessarily given 
by those in eq.(j4]). While the mass matrix in the form of eq.([^ provide sufficient condition for 
5 = ± 7 r /2 and ^23 = 7 r /4 when S 13 and sin 5 are not zero. The simplicity of the above mass matrix 
may serve as a good starting point to understand possible underlying theory. In fact it has been 
shown that the above neutrino mass matrix is a consequence of imposing a symmetry of the form 
e —)■ e, /i and r exchange with CP conjugation discussed by Grimus and Lavoura in Ref. [^, which 
we will refer to as the Grimus-Lavoura symmetry (GLS). 

In this work, we study realizations of <5 = ± 7 r /2 and ^23 = 7 r /4 in type-II seesaw model with 
A 4 , flavor symmetry. Models based on ^4 symmetry has been shown to be able to provide a good 
scenario to achieve this[ll|, 12|. In ^4 models, the charged lepton mass matrix Mi is diagonalized 


from left (rotation on left-handed charged leptons) by the characteristic matrix Ui for ^4 symmetry 
model buildings [^, 


M 

Ml = UimiUr , G; = ^ 1 a; , (8) 

V 3 1 u; ) 

where u = exp{i27i/3) and = exp{iA7r/3). Ur is a unitary matrix, but does not play a role in 
determining Vpmns- We will not specify its form here. 

If neutrinos are Majorana particles, the most general mass matrix is of the form 

( Wi X y \ 

X W2 z \ , (9) 

y z W 3 J 

which can be diagonalized by unitary matrix W, M^, = V^rh^Vj. 
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The mixing matrix Vpmns is given by 

VpMNS = U^Vy . ( 10 ) 

In the basis where charged lepton is diagonalized, the neutrino mass matrix is of the form given by 
eq.(0 


a = -{wi + W2 + W3 + 2 {x + y + z)) , 

d + i'y = ^{wi + 0JW2 + + 2{<jj'^x + ojy + z)) , 

o 

d — i'y = -{wi + u‘^W2 + 0JW3 + 2{bjx + oj'^y + z)) , (11) 

3 

c + ip = - {wi + uj‘^W2 + 0JW3 — ujx — oj'^y — z) , 

3 

c — ip = -{wi + UW2 + — u'^x — (jjy — z) , 

3 

b=^{wi + W2 + W3-{x + y + z)). 


If one imposes the GLS on the neutrino mass matrix, all parameters in the set P = {wi,x,y,z) 
are dictated to be real, and will predict (ill, fl^ 6 = ±7r/2 and 623 = vr/4. Therefore in A 4 model 
building for neutrino masses with 6 = ±7r/2 and 623 = vr/d, it is essentially to make sure that Ui 
is of the form given by eq.([H]) and require the resulting mass matrix to satisfy GLS. Note that in 
this case since the parameters in set P are all real, the complexity of the mixing matrix is purely 
due to the appearance of u and When the GLS is broken by allowing the parameters in P can 
be complex, there are more sources for GP violation and the model does not predict 6 = ±7r/2 
and 023 = 7r/4 automatically. This points a way to modify the predictions to £t data should future 
experiments will find 6 and 623 to be deviate signihcantly from —7r/2 and 7r/4. We will study both 
cases with the parameters in set P to be real and complex in the rest of the paper. 


II. TYPE-II SEESAW MODEL WITH A4 SYMMETRY 

We now construct a Type-II seesaw model [l3| with A4 family symmetry to realize the forms 
of mass matrices in eqs.([6]) and ([7]). A different model based on Type-II seesaw with A 4 has 
been constructed to realized tribi-maximal neutrino mixing[l^. In our model, the left-handed 
lepton doublet Ip and right-handed charged lepton singlet Ip have the following standard SU{3)c x 
SU{2)i X U{1)y gauge and A 4 family symmetry properties 

h : (1, 2, -1)(3) , Ir : (1,1, -2)(1 + l" + T) , (12) 

where the first three numbers in the first bracket indicate the SU{3)c, SU{ 2 )r and U{1)y trans¬ 
formation properties. The numbers in the second bracket indicate the A4 representations. 

To obtain desired mixing pattern, the Higgs sector is enlarged to have two types of Higgs 
doublets, p and and two triplets, A and y for neutrino masses. They transform as 

P : (1, 2, -1)(1) , $ : (1, 2, -1)(3) , A°': (1, 3, -2)(1 + T + l") , y : (1, 3, -2)(3) . (13) 
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The Lagrangian responsible for the lepton mass matrix is 


L = yelL^lR + y^,lL^l\ + yrlL^l]i 

+ + Ylh^ll + yX^'II + vJlxII + H.C. 


( 14 ) 


If the structure of the vacuum expectation value (vev) is of the form < >= uf, < Xi >= 

>= v^’ , one obtains the charged lepton and neutrino mass matrices Mi 


< (j) >= ^ 0 , and < A°’ 
and My as 



I/e 0 0 

0 0 
0 0 yr 


, My = 


Wi X 
X W2 


y 


z 


(15) 


where 


wi = Yy\ + yIv^ + Y”v'^ , 

W2 = Yy^ + u:Xv^+^y'^vl, 

w, = fx + , 

X = yuv^ , y = yX , 2 ; = yX 


(16) 


If the vevs of < < 1 ) 12,3 >= '^f 2 3 equal to v^, the vev structure of the Higgs helds breaks 

H 4 , but left with a residual symmetry generated by {I,c,a}. Here c, a are H 4 group elements 
dehned in Ref. 1^. This will lead to the charged lepton mass matrix of the form 


me 0 0 

Mi = Ui \ 0 m^ 0 
0 0 m.. 


(17) 


where ^ \/3yr,e,rr'^. 

As long as Z 3 residual symmetry in the model is not broken, we have one of the key ingredients 
in realizing the form of charged lepton mass matrix in eq.([ 8 ]). In the basis where the charged lepton 
mass matrix is diagonalized, the neutrino mass matrix will be the same form in eq.(| 6 ]). 

In the rest of this paper, we will study consequences related to the above mass matrix. We will 
study conditions on the model parameters imposed by GLS, and phenomenologically acceptable 
models can result in this class of models. 

Before carrying out detailed analysis for neutrino mixing, we would like to point out that this 
model can easily accommodate data in the quark sector if one assigns the left- and right- handed 
quarks Ql, Ur ad Dr as A4 singlet “1”. In this case the Yukawa couplings for quark masses and 
their mixing are given by 


L = QlYucPUr + Wl^r + H.C. . (18) 

This will give, in general, arbitrary 3x3 up- and down- quark mass matrices after cj) develops 
a non-zero vev with Mu,a = Yu,dV(f). These matrices will have no predictive power for quark 
masses and their mixing, but have no problem in accommodating experimental data. 
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III. MASS MATRIX PREDICTING <5 = ± 7 r /2 AND 623 = 7 r /4 


If the parameters in the set P are all real, the resulting mass matrix is of the form given in 
eq. Therefore this model predicts 


5 


±- , O 23 
2 


71 

4 ■ 


(19) 


The above prediction can also be obtained by studying the mixing matrix Vpmns is eq. ffTOj) . For 
an arbitrary complex symmetric matrix M^, the matrix I 4 diagonalizes it, has the most general 
form 


K = dtag{l, ^” 13 , O 23 , S) ■ 


( 20 ) 


Here we have absorbed a Majorana phase matrix diag{l, on the right of the above equation 

into the eigen-masses 

If the parameter set in P are all real, the mass matrix can be diagonalized by an orthogonal 
mixing matrix Vy = 6 * 13 , ^ 23 , <^ = 0), and are real. Depending on conventions, here Sij and 

Cij need not to be normalized all to be positive. We have the Vpmns given as the following 


— S12(C23 —S23)+C12(C13 —^13(c23+.g23)) C12(c23 —^23)+^12 (ci3 —^13(c23+S23)) gl3+C13(c23+^23) 

V3 _ __ _Vf_ _ 

—S12iPc23—‘^S23)+C12ici3—Sl3{u]C23+Ps23)) Cl2(Pc23—i^S23)+Sl2ici3—S13(u}C23+Ps23)) S13+C13(lt^C23+1^^^23) 

Vs __ Vs 

-Si2{LdC23-Vs23)+Ci2{ci3-Sl3{Vc23+LdS23)) Ci2{u}C23-V S23)+Sl2ici3-Sl3{V C 23 +LJS 23 )) Si3+Ci3(Vc23+UJS23) 


\ 


( 21 ) 


V3 


F3 


F3 


The above matrix has the property that, \VpMNS, 2 i\ = |HpMAfS, 3 i|- Rewriting the above into the 
standard form of the mixing matrix, this leads to[9| 

923 = — , si3Cos(5 = 0. (22) 

Since experimentally S 13 7 ^ 0, then cos <5 = 0 which leads to 5 = ± 7 r/ 2 . 

The other angles can be determined by 

513 = IR 13 I , tl 2 = |tT-7 , (23) 

I *^111 


where tij = Sij/cij. 

The condition whether 6 takes -l- 7 r /2 or — 7 r /2 is determined by using the Jarlskog parameter (l^ 
J = Ci3Si2Ci2S23C23Si3sin5 = /m(l/ii 1 / 12 ^ 2 *^ 22 ) to obtain 


(l-|W3p)Jm(RnRi2^*i^22) 

IR11IIR12IIR23IIR33IIR13I 


(24) 


Note that if 5 7 ^ 0 and/or are not real, in general |VpMAfs, 2 i| is not the same as |VpMvs, 3 i|- ^23 
and 5 will not necessarily be 7 r /4 and ± 7 r/ 2 . One should be cautious about this. This uncertainty 
actually can be turned into a good use to provide a possible way to modify the predictions for 623 
and 5. The value for 6*23 is determined by 
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The angles 613 and 612 are determined by eq. fl2^ . and the CP violating phase 6 is given by eq. fl24|) . 

The above discussions provides a way of looking at how the reality of the parameters in the set 
P leads to 5 = ±7r/2 and 623 = vr/4. It is essential to have the parameters in the set P to be all 
real. The complexity of the parameters can appear in the Yukawa couplings, in the vevs, and also 
in places where cu* appear in neutrino mass matrix. One needs to see with what conditions they 
can be made real. 

To make the Yukawa couplings and scalar vevs real, one can require the model Lagrangian to 
satisfy a generalized CP symmetry under which 


{le,L 

(A^ 


A 


A") ^ (A°t , A't , A"t) , 


(xi , X2 , X3) ^ (xl , X3 ) xl ), 


(26) 


and all other helds transform the same as those under the usual CP symmetry. Here the superscript 
CP in the above indicate that the helds are the usual CP transformed helds. 

The above transformation properties will transform relevant terms into their complex conjugate 
ones. Requiring the Lagrangian to be invariant under the above transformation dictates the Yukawa 
couplings to be real. The same requirement will dictates the scalar potential to forbid spontaneous 
CP violation and vevs to be real. One, however, notices that the parameters tC 2,3 in eq. flTB]) are 
in general complex even if the Yukawa couplings and the vevs of the scalar helds are made real 
because the appearance of cu*. To make them real to satisfy GLS, it is therefore required that 

+ ujY”v'^) = Im{ujY^v'^ + u/y”v'^) = 0 . (27) 

The above can be achieved by the absent of the scalar helds A'’” in the theory or YJua' = Yu'^A"- 
We will consider examples for each case in the next section. 

From the above discussions, we know that in general the mass matrices obtained can accom¬ 
modate values diherent than vr/d and ±7r/2 for 623 and 6 respectively, we will not carry out a full 
numerical search analysis for parameter spaces, but to take some simple cases to show how the 
modihcations arise when GLS is broken and how future experiments can test this class of models. 
These models, due to additional constraints, will have some additional predictions than the general 
one. 


A. Model A: Z2 residual symmetry for neutrino mass matrix 


If vev of X 2 component of y is non-zero, but the vevs of yi ^3 are zero, the vev structure breaks 


A 4 down to a Z 2 generated by {l,r 2 }. Here r 2 is an A 4 group element dehned in Ref. 
charge lepton and neutrino mass matrices are given by 


15 . The 


/ me 

0 

0 ) 

1 j 

/ Wi 

0 

y \ 

0 

m^ 

0 

, M,= \ 

0 

W2 

0 

VO 

0 

rur ) 

1 1 


0 

Wz / 


(28) 


The parameters in the set Wi, and y are in general complex. 

If the helds. A' and A" are absent from the model, one would have Wi = W 2 = W 3 = Y^v^. The 
resulting mass matrices lead to the well known tribi-maximal mixing 17| which had been the focus 
for A 4 symmetry model buildings before 15l. Il8l-l20|. In this case, 643 = 0 and phase 6 is zero. Since 
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experimentally S 13 has been determined to be non-zero, this model is ruled out. Modihcations have 
to be implemented j2ll-l^ to accommodate data. A way out is to keep in the model 

as discussed in the previous section. 


We now summarize the main results follow Ref. 11 for this model. This also services as an 


outline how analysis can be carried out for relevant models. Letting M^, be diagonalized by 14, one 
has 


= V^rhJ/J . (29) 

14 in this case can always be written in the following way with Majorana phases to be absorbed in 
to the neutrino masses m*. 


14 = 1414 , with 14 



(30) 


where I 4 is a diagonal matrix diag{l, 1, e*^) with tanp = Im{yw\+y*w:i) jRe{iiw\+y*w^). s = sin6* 
and c = cos 6*. 

Expressing the mixing angle 9 in terms of the model parameters, we obtain 


tan 20 


2\ywl + w:iy* 

kip - ksP 


(31) 


Using rhy = V}V^Mi,V*V* derived from eq. (|2^ . we obtain the Majorana phases Oj of m* in 
terms of the model parameters as 

ai ,3 = Arg{wi{l ± cos 20) -|- =F cos 20) ± 2 sin 20|/e“*'’ , 0.2 = Arg{w 2 ) ■ (32) 

One can always normalize the Majorana phase a 2 to be zero without of generality. In this basis, 

the phase of W 2 is also zero. Although the Majorana phases can be expressed in terms of the model 
parameters, since there are more number of parameters than the mixing angles and eigen-masses, 
there is no prediction for the Majorana phases. We will not discuss them in our later numerical 
analysis any more. 

The mixing matrix Vpmns can be, in general, written as 

. / c + 1 cRP — s \ 

VpMNS = UlVpVy =\ c + ojsRP u? ojcRP - s 1 , (33) 

y c -I- oj'^sRP u ~ ^ ) 


Using eqs. (l2^ . (12T)) and ([25]), we find that 


S 12 = 


\/2(l + 


cs cos 




523 — 


(1 -f- CS COS p -|- \/3cs sin 
\/2(l -1- cs cos p)^ 


( 1 - 2 cscosp)i /2 

-. (34) 


and 


sin (5 = (1 -|- 


4cVsi^ 1 / 2/3 
(c2 - S42 ^ ^ 


3c^s^sin^p 1/2 f -1 ) if > s^ , 

(1 -f cs cos py I ’ if s^ > . 
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(35) 










From the above, one clearly sees that if sinp is not zero, |(5| and 623 deviate from 7 r /2 and 7 r/ 4 , 
respectively. In the limit p goes to zero, however, the above recovers the results with real parameter 
set P with 6 = ± 7 r /2 and O 23 — 7r/4. 

There are two interesting features for this model worth mentioning. One of is that \Ve 2 \ to be 
l/\/3 which agree with date. S 12 is always larger or equal to l/\/3 which is a decisive test for this 
model. Another is that although the Dirac phase 5 depends on the phase p, the Jarlskog parameter 
J which is independent of p given by J = — (c^ — s^)/6\/3. This implies that CP violation related 
to neutrino oscillation is still purely due to intrinsic CP violation. This model can be made in 
agreement with data at la level. 

If p = 0 and c = s = l/\/2, the mixing pattern is the tribi-maximal. However, if p is not zero, 
even if c = s = Si 3 can be non-zero,S 12 and S 23 are also modified from their tribi-maximal 

values 


1 1 \/3sinp (1 —cosp)^/^ 

^ (2 + cosp)V2 ’ ^23 = ^(1 + 2 + cosp^ ’ ^ 7! 

J is exactly zero which implies sin 5 = 0. 


(36) 


B. Model B: mass matrix with tci 2,3 = 0 


Without 5 Wi^ 2,3 = 0. In this case all the three vevs of y should not be equal in order to fit 
data on the neutrino mass-squared differences = |m 2 p — |n7.ip and Arug^^ = jmsp — |mip. 

The Z 2 residual symmetry mentioned previously is also broken in this case. If CP is spontaneously 
broken in the Higgs potential, nf may be complex too. Therefore the parameters in the set P : 
{x,p,z} = {|x|e*°'“’, |p|e*°''', are in general complex. One can always choose, without loss of 

generality, to rewrite the neutrino mass matrix as = VaMyVj 

/O |x| |p|\ /I 0 0 \ 

= |x| 0 1^1 , 14 = 0 0 , (37) 

\ \y\ \z\ 0 J \0 0 J 

where ai = (a^ — — o'y)/2, (T 2 = (Tz — Cy and a 3 = a^, — ay. ai can be absorbed into lepton helds. 

We will neglect it in the following analysis. 

Let us discuss now how the model parameters can be determined by data. First of all since M^, 
is a real matrix, the eigen-values m* are all real and fhi + m 2 -|- m 3 = 0. The absolute value of 
neutrino mass |mj| = \fhi\ can be solved as functions of Am|]^ and Am|^. For the values allowed 
by experimental data, the solutions are given by up to a overall sign 

(Am^i — AA^)7^^2 i ~ + 2AA2 

\/3(Am2i + Am|]^) 

(Amgg -|- AA^)7^^2 i ~ + 2AA2 

\/3(Am2i -l- Am|]^) 

AA^ = \J{Arnhy + Am^j^Arngg -|- (Am|g)2 (38) 

One can easily check that multiplying a ” to the above one obtains another solution. The overall 
sign can be absorbed into lepton field phases redefinition. We will use the above normalization for 
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the overall sign. In the numerical study later, we will show that only inverted hierarchy is allowed. 
Therefore, mi ,3 > 0 and m 2 < 0. Since the eigen-masses are all real, there is no Majorana phases 
in this model. 

Since is real, it can be diagonalized by a orthogonal matrix. Assuming U (^13 , ^13, ^23, (5 — 0 ) 
diagonalizes M^, one can express the fact that the “ 11 ”’ and “ 22 ” entries are zero as 

Ci 3 (miCi 2 + hisSia) + ^3513 = 0 , 


hfi(si 3 — ii2i23){ii2 + i23Sl3) — m2(tl2Sl3 -|- t^23)(l ~ ^12^^23513) + '^3^i23 — 0 , (39) 

Cl2 

where iij = Sij/cij. 

From the above two equations, one can express S 12 and §23 as functions of §13 and m^. Since 
rhi are known, the elements in Vpmns = UjVfjVi, can be expressed as functions of (T 2 , (T 3 and § 13 . 
Using the relations between the elements Vij of Vpmns with Sij and S in eqs. fl2^ and fl2T|) . one can 
vary (T 2 , (T 3 and §13 to see if the resulting values for 512 , 13,23 are in the allowed values, and obtain S 
by ea.fED. 

Similar analysis had been used to rule outjl^ the simple version of Zee model where the mass 
matrix for neutrinos is M^, in the diagonalized basis of the charged leptons. Here the appearance of 
Ui may save the model. We hnd that, unfortunately, that with real parameters in P which implies 
( 72,3 to be zero, 6 and 623 are predicted to be — 7 r /2 and vr/d, but there is no solutions within the 
allowed ranges which can predict S 12 to be consistent with data. With complex parameters in P, 
the values for 512 , 13,23 can be in agreement with data, it is, however, not possible to get 6 to be 
close to — 7 r/ 2 . More numerical details will be provided in the next section. 


IV. NUMERAL ANALYSIS 


In this section, we compare experimental data with the model predictions for the mixing angles 
and CP violating phase. There are several global hts of neutrino dataj^, |^. The latest £t gives the 
central values, la errors and the 2 a ranges as the following!^ 


S/ti 

NH 

2 cr region 0.0 ~ 2.0 

IH 1.48 ±0.31 

2cr region 0.00 ~ 0.09&0.86 ~ 


*12 

0.323 ±0.016 
0.292 ~ 0.357 
0.323 ±0.016 
2.0 0.292 ~ 0.357 


4/10-^ 
2.26 ± 0.12 
2.02 ~ 2.50 
2.29 ±0.12 
2.05 ~ 2.52 


523 

0 e;fi7+0-032 
U.OO^_Q 324 

0.414 ~ 0.623 

n e;7Q+0-025 
U.O(0_o.o3g 

0.435 ~ 0.621 


and the corresponding values for mass-squared differences are given by 


NH 

2 a region 
IH 

2 a region 


Aml[W-^]eV^ 

7.60;»;!» 

7.26 ~ 7.99 

' -cu.o.is 
7.26 ~ 7.99 


|Amli|[10-3eW 
2 48+°-°5 

^•^O-0.07 

2.35 ~ 2.59 


+ 0.05 

- 0.06 


2.38 
2.26 ~ 2.48 


(40) 


(41) 


Here NH and IH indicate neutrino mass hierarchy patterns of normal hierarchy and inverted 
hierarchy, respectively. We will use the above data for comparison. 
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A. Some general comments 


The model discussed achieved a natural way to have the CP violating angle 6 to be ± 7 r /2 
and 623 to be 7 r/ 4 . The condition is to have all the parameters in the set P to be real. This 
requirement can be taken as a consequence of having CP violation only caused by intrinsic source 
of CP violation. The condition for choosing 6 = — 7 r /2 can be determined using eq. fl24p which 
requires /m(l/iiP,* 2 m^ 22 ) to be negative. 

The value — 7 r /2 predicted in the model is in agreement with IH within la range. Although for 
NH case 6 is outside of la range, there is no problem with 2a range. For S 23 , the model predicts 
S 23 = 0.5. This value is outside of la range for both the NH and IH cases. However, they are, 
again, in agreement with data within 2a. Other parameters, Am?^, S 13 and S 12 can be easily htted 
to within la level of current experimental data. 

With complex parameters in P, 6 = —tt/2 and ^23 = 7 r /4 can also deviate away from —ti/2 and 
7 r /4 with sin5 determined by eq. (l24|) and 623 determined by eq. fl25l) . One can hnd solutions where 
^ 12 , 13,23 and 5 to take the central values from current data. 

As the general model is expected to be able to ht data, it may be more instructive to analyze 
some simplihed versions than just providing with numbers. We provide more details of Model A 
and Model B discussed earlier next to see how additional assumptions restrict the level of model 
agreement with data. 


B. Model A predictions 

With real values for Wi and y in Model A, by adjusting the values, tCj, and y both NH and IH 
mass patterns can be obtained. The predictions for 5 and 623 are ± 7 r /2 and 6*23 = 7 r/ 4 , just like the 
general case discussed in the previous section. Since 5 should be close to — 7 r/ 2 , one should take 
the parameter space so that (? > s^. In the model S 13 = (1 — 2cs)^/^/\/3 is not predicted. But one 
can use information from S 13 to £x cs = 0.497 ± 0.018 to predict s ^2 = 0.334 ± 0.004 for both NH 
and IH cases. This is in agreement with data within la. Note that = ( 512013 )^ = 1/3. It agrees 
with data within la. 

It is remarkable that neutrino mixing matrix in this model with just one free parameter can be 
in reasonable agreement with data. This may be a hint that it is the form for mixing matrix, at 
least as the lowest order approximation, that a underlying theory is producing. One should take 
this mass matrix seriously in theoretical model buildings. 

If the parameters in the set P are complex, and therefore a new phase p appears in the model. In 
this case, the new parameter p can be used to improve agreement of the model with data. In both 
NH and IH cases, 5 and S 23 can be brought into agreement with data at la level. As an example, 
we take the largest value of cs so that S 13 takes its lower la allowed value, and then varying cos p 
to obtain the upper la allowed value. This hxes cs and cosp to be 0.468 and 0.992, respectively. 
With these values, S 23 and 5 are determined to: 0.534 and 1 . 4267 r, respectively. These values are 
in agreement with data at la level. 

For the case with c = s, the model is more restrictive. In this case sin 5 = 0. With more precise 
data on the CP violating angle 6 , this may rule out this simple case with high conhdence level. If 
p = 0, the model is already ruled out at high precision from S 13 measurement. However, with a 
non-zero p, the mixing angles can still be made to in agreement at 2a level. In Fig. [T]we show S 12 , 
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S 23 , 5 i 3 as functions of p. When chose cosp = 0.93, we can get S 12 = 0.58, S 23 = 0.78, S 13 = 0.15 
which agree with the experimental data within 2 (T range. 





FIG. 1: S12, S23, S13 as functions of cosp with c = s = V^l2 compared with allowed ranges for 
(•Sis, 'S23, S12) at Icj and 2cj, respectively, given by (0.147 ~ 0.155, 0.731 ~ 0.773, 0.554 ~ 0.582) 
and (0.143 ~ 0.159, 0.660 ~ 0.788, 0.540 ~ 0.597). 

More precise experimental data are required to distinguish the model with complex model pa¬ 
rameters from that with the real parameters and other models, or to rule out the above simples 
completely. 


C. Model B predictions 

With real parameters in P in this model, the predictions, 6 and 623 are the same as the general 
case. Using eq. fl39l) . one can choose §13 to be the parameter to fit all data. We find that no solutions 
can simultaneously bring S 13 and S 12 to be in 2a allowed region compared with data. In Fig. [2] 
we show S 13 , S 12 as functions of § 13 . Here, when S 13 agree with the data in 2 cr range, the allowed 
region for S 13 are (—0.594 ~ —0.565)&(—0.012 ~ —0.010)&(0.010 ~ 0.011), but when S 12 agrees 
with the data in the 2a ranges, the allowed region for §13 are only (0.806 ~ 0.774), they have no 
overlap region for § 13 . This model is therefore ruled out. 




■^13 


•^13 


FIG. 2 : S13, S12 as functions of S13. 

With complex parameters in P, there are solutions of S 13 and CTj so that to make S 23 , S 13 and S 12 
to be consistent with data at la level. However this requires one of | sinujl to be away from zero 
which alter the prediction for 6 signihcantly away from —7r/2, although 623 is still close to 7r/4. In 
Fig. 131 we scan the parameter space —1 < S 13 < 1, 0 < (T 2 < 27r, 0 < as < 2%, and obtained the 
corresponding values of S 12 , S 23 , S 13 that agree with data within la and 2 a region respectively. And 
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in la range, the CP phase will be constrained in sin 6 ~ (—0.5, 0.5), but in 2a region, CP phase is 
free and can be in the range of 0 to 27r. 


2cr range 



0.154 


FIG. 3: The scanning points of the three free parameters in la and 2a region respectively. 

For example with sis = 0.57, a 2 = crs = 267 r/ 25 , one can obtain that sis = 0.15, S 23 = 0.74, 
S 12 = 0.57. These mixing angles are within the Icr allowed region, but it will have a CP phase 
(5 = I.OOtt signihcantly away from 37 r/ 2 . A precise determination of the CP violating phase is 
needed to rule out the above model. 


V. CONCLUSIONS 

We have constructed theoretical models for neutrino mass matrix in Type-II seesaw with A 4 
family symmetry. The models we constructed naturally predict that the CP violating phase 6 is 
equal to ± 7 r /2 and at the same time the mixing angle 023 is 7 r/ 4 . The reality of the parameters 
can be achieved by imposing a generalized CP symmetry in the Grimus-Lavoura symmetry limit. 
When the generalized CP symmetry is explicitly broken, the Yukawa couplings can be complex, the 
model predictions for 6 and 6*23 can be signihcantly modihed. Two simple scenarios. Model A and 
Model B, are analyzed in detail to show how the modihcations arise and how future experimental 
data can test this class of models. 

The Model A has a characteristic prediction that |Ue 2 | = l/\/3. This model can accommodate 
experimental data at la level. It can be taken as the lowest order neutrino mass matrix for future 
theoretical model buildings. When Yukawa couplings are complex, the CP violating phase 6 and 
the mixing angle 623 can be away from — 7 r /2 and 7 r/ 4 . The crucial test for this model is to measure 
whether \Ve 2 \ = l/\/3 holds to high precision. 

In Model B, the diagonal entries of neutrino mass matrix are all zero in the weak interaction 
basis. This implies that the neutrino masses can be determined by the known neutrino mass-squared 
differences. We hnd that this model can only accommodate inverted neutrino mass hierarchy. The 
mixing angles S 12 , S 13 and S 23 cannot simultaneously be in agreement with data at 2 a level if the 
Yukawa couplings are all real. With complex Yukawa couplings, the mixing angles can be brought 
into agreement with data at la level, but the CP violating angle <5 will be signihcantly away from 
— 7 r/ 2 . This provides a crucial test for this model. 

At present, experimental data on neutrino mixing seem to hint that 6 and 623 to be close to 
—tt/2 and 7 r/ 4 . Theoretical models which can naturally achieve such predictions are interesting to 
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study. The models we have constructed have many novel properties and can be tested. We will 
have to wait future experimental data to tell us more whether this class of neutrino mass matrix 
will survive. 
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